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Using a simple analytic expression for qq, gg — > tt — > bW + bW~ — » blvibl'vli with 
the interference effects due to the polarizations of the t and t, we demonstrate how 
the effects can be measured at Tevatron at 3cr level. 



With the upgrade of Tevatron and the improvements of its two collider detectors, the data related 
to the production of the top quark are going to accumulate very quickly. Since top quark decays very 
, rapidly after its production, it is expected that the spin information of the top quarks is preserved in the 

decay process. For lighter quark, this spin information is always smeared by the hadronization effect. 
Therefore, if such effect can be observed it will be the first time that one can observe the spin of a 
bare quark directly. However, it is challenging to investigate how one can determine the various detail 
properties of top quark in the complex hadronic collider environment, 
■^-j- ' Here we report our investigate an interesting physical consequence of the fact that the top quarks 

(S) | that are produced and decayed are supposed to be spin 1/2 particlecla. One of the important effect of 

the polarizations of unstable particles are that the different polarized intermediate states can interfere. 
In this sense, the observation of interference effects provides a unique possibility of direct observation of 
the spin of a quark. 
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[ We use a simple analytic resulltl for the differential cross section of qq and gg — > tt — ► bW + bW~ 

blvjbl'vy based on an analytic helicity technique developed in ref. a. The decays of W bosons and top 
quarks are taken into account in the narrow width approximation. The contribution due to off-shell top 
quarks or offjshell W bosons are negligible. The interference effects discussed here was also considered 
before in ref. u, however, it was studied only numerically and was done in a rougher approximation. The 
analytic expressions obtained here can also be easily adapted to the leptonic colliderfl. We shall also 
demonstrate that such interference effect can be detected at 3<r level at Tevatron with Main Injector if 
one uses proper observablefl. 
Oh The polarization density matrix for the process can be split into two main sections and written as 
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Here the initial the final polarizations are summed over. N s = A(N 2 — l) 2 for gg initial states; and N s — 
AN 2 for qq initial states. Ii q {ri) = —i(rf — m 2 ,)~ 1 represents the polarization independent components 
of the top quark propagators. The helicity informations are included in the remaining density matrices. 
When no subscript is included for a given particle with spin in P, it implies that its helicity is summed 
over. The polarization density matrix P{i\ 1 \ 1 i\ 2 \ 2 — > ^a ! a^a 4 aj 1 ) represents the production of the top 
quark pair via the processes qq — * tt or gg — > tt. P(t\ 3 x ?a 4 a^ ~ > bbllvv) is the polarization density 
matrix for the decay of the top quark pair into b quarks and leptons. This density matrix can be further 
splitted into the product of the decay density matrix of tt into W + W~ bosons and bb pairs, and the 
decay density matrix of W + W~ boson pair into a pair of fermions each. 

The polarized density matrix of the process t\ 3 y ^a 4 a^ ~~ * bbllvv can be written as 

p ( *a 3 a' 3 W 4 - m+i-vv) = p ^y 3 ^ w tx>b)P(h iK ^w- 6K b) 

x \n w ( Pl )\ 2 \n w ( P2 )\ 2 P(W+ - i+v)p(w- 6K - i-v). (2) 
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Simple analytic expressions for all^he polarization density matrices appeared in the first hand side of 
above equation can be found in Reftl. These expressions can be combined together with the narrow width 
approximation for the W boson propagators and obtain 
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The helicities of final state fermions are summed over. We shall use four vectors fci, U and qi to 
represent the momenta of the neutrinos, the charged leptons and the b quarks respectively. The subscript 
1 is used to label the momenta related to the decay products of the top quark, while the subscript 2 is 
used for those related to the decay of the anti-top quark. e lM (A5) represents the polarization vectors of 
W + with helicity A 5 , while e 2 u{\6) represents those of the W~ with helicity A 6 . In our notation, the 
information about the top helicity is carried by the Pauli matrix m and cto(= 1) with component (A3 A3). 
For the anti-top, we used the conjugate matrix a a — a 2 a a a 2 with component (A4A4). The tensor e a (Yi), 
(a = 0, 1, 2, 3), were introduced in rem to evaluate the product of fermion wave functions, u(r, X)u(r, A'), 
in the amplitude squared for the cross section. 

With this, one can then proceed to glue together the density matrix for the production, which has 
been given in RefH, and those for the decays of the top quarks. For the top and anti-top propagators, we 
shall use also the narrow width approximation. These provide another two delta functions. After that, 
we can integrate out all the delta functions in the phase space for 6 final state particles and obtain the 
differential cross section 
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where 



+ i 2 , with i\ and i 2 defined as the four- momenta of the initial partons and d 2 (a — > be) is 
the phase space for a particle with four-momentum a decaying into particles with four-momenta b and c. 
N is the number of color and the indices ii designate the initial particles. For the gg — > it case, we have 
the result 
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and for the qq 
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The s, t and u are the Mandelstam variable, i? M 
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s/4m 2 and x = s 2 /2(t — m 2 ){u — m 2 



To show the spin effect of the quark, we shall compare with the similar calculation in which top 
quark helicities are averaged over both in the production density matrix and in the decay density matrix 
independently before we join them together. We shall refer to this case as the "spinless case" . In the usual 
Monte Carlo simulation for the experiments, this is indeed the approximation taken by the program. The 
cross section for the spinless case can be obtained by substituting I gg and I qq with 



jafi 
99 



and 



1 



N 2 -l 



1 '2 



1 

X 



2r 



(8) 



m 



N 



2N 



1 II 

x 



2 7 2 



(9) 



For the last two equation, we should emphasize that the cross section without the spin correlation for 
the top quark is the one that has been used in the literature and by most of the experimental event 
simulators. 

To measure interference effect, we need to find observables that correlate the kinematic variables as- 
sociated with the particles from the top quark decay and those from the anti-top quark decay. Intuitively, 
it is not obvious what is the best correlated observable which can probe this interference effect. For the 
rest of this paper, we shall present our attempts in this direction. 

To probe the interference effect, we shall look for observables related to the final decay products, 
instead of using directly the top momenta which need to be reconstructed. The easiest observables to 
try are those related to the two charged leptons (electron or muon) from the leptonic decays of two W's. 
The advantage is that these particles can be easily observed in a hadron collider. On the other hand, the 
dilepton decays has lower branching ratios, which reduces the statistics. For this case we had investigated 
a the effect of interference terms on the distribution of (1) their total energy Ei t + Ei 2 , (2) their total 
z-momentum l\ z + l<z z , (3) an orthogonal combination l\ z — l^z, (4) the cosine of the angle between these 
two charged leptons and (5) the asymmetries Ap to be defined later. 

Here we shall simply report the observables that we think are most promising. We shall demonstrate 
that it is possible to observe such asymmetries at 3u level with Tevatron II. The lepton correlated 
asymmetry, Ap, with respect to a plane P passing through the interaction point is defined as follows. 
Let l\ be the momentum of the charged lepton associated with top decay evaluated in the top rest frame. 
Similarly, let I2 be that associated with anti-top decay evaluated also in its own rest frame. Then, define 
Ap = (N + — N-)/(N+ + N-) where N + is the number of events in which both li and I2 lies on the 
same side of P while iV_ is the number of events with both l\ and I2 lying on the opposite side of P. 
By choosing different P, one can construct different asymmetries. These asymmetries vanish when the 
effects of both top spin and W spin are ignored (or averaged over) and they remain small even when the 
W spin effects are included. This property makes them the ideal candidates for the observation of the tt 
spin correlation effects. 

In order to measure the lepton correlated asymmetries it is essential to make full reconstruction 
of the top dilepton events. We assume that the top mass will be well measured in Tevatron Run I 
and Run II. For dilepton candidate events, we further assume that the missing transverse momentum 
measured is equal to the sum of the missing transverse momenta of the two neutrinos associated with the 
dilepton. Contribution to the missing transverse momentum from other neutrinos in the event reduces 
the efficiency of the reconstruction and lowers the signal-to- noise ratio but does not spoil the observability 
of the asymmetries we consider here. 

The analysis of the productions and the .decays of tt in hadronic collider with and without spin 
correlation effect can be found in the literaturolj with various degrees of sophiscations. We shall use t 
simple analytic differential cross sections for qq and gg — + tt — > bW + bW~ — ► blvibl'vi' provided in Re 
Using these formulas, we simulated the tt production and decay employing the event generator PYTHIA 
5.7. The algorithm of our reconstruction of the top dilepton events has been described in detail in Ref. □. 
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In principle, one can also look for the W correlated asymmetries for the two W's from tt decays in 
the reconstructed top dilepton, lepton-jet and jet-jet events. The predicted W asymmetries can be easily 
obtained from Refu. They are no more than a few per cent and about an order of magnitude smaller than 
some of the lepton correlated asymmetries. The enhancement of the lepton correlated asymmetries over 
the W correlated asymmetries] is one feature that favors the observation of top spin correlation effects 
through the dilepton channel. 

In hadron colliders, tt are produced either by quark-antiquark (qq) annihilation or by gluon-gluon 
(gg) fusion. The lepton correlated asymmetries for these two cases typically have opposite sign and with 
different magnitude!]. At Tevatron energy, the number of qq produced top events are roughly eight times 
that of the gg produced events. This ratio decreases with energy and, at LHC energy, the number of 
gg produced top events are roughly four times that of the qq produced events. For the asymmetries we 
considered, we found that the substantial asymmetries that can be observed at Tevatron become much 
smaller at the expected LHC collider energy. This is the result of accidental cancellation between the 
contribution of the qq production channel and that of the gg production channel. We have checked that 
the contribution due to the qq channel alone is indeed not so small. This cancellation could be a generic 
feature which may make LHC unfavorable machine to look for top spin correlation effect. For the same 
reason, the next linear colliders should be an ideal machine for observing such correlation due to the 
absence of such cancellation. 

As examples, we shall discuss the asymmetries with respect to the following planes defined in the 
tt center of mass frame: (1). tt production plane (defines asymmetry A\)\ (2). the plane perpendicular 
to the production plane and contains the top (asymmetry A 2 ); (3). the plane perpendicular to the two 
previous planes (asymmetry A3); (4). the plane normal to the beam direction (asymmetry A 4 ). These 
planes are chosen as samples to demonstrate the possibilities. The optimal choice will depend on the 
details of the detectors and clearly need further studp. 

The typical trigger for top dilepton events, namely, pt > 20GeV for leptons, pr > lOGeV for 6-jet, 
missing transverse energy > 25GeV and rapidity \r]\ < 2, was applied. Afterwards, the reconstruction 
algorithm we described earlier was carried out. A typical hadron calorimeter energy resolution of 70%/ \[E 
was used as b - jet energy smearing and the missing transverse energy was Gaussian smeared with a 15% 
standard deviation. The effect of including contribution of other neutrinos in an event to the missing 
transverse energy was investigated! To isolate the effect of the bias originated from reconstruction 
algorithm to the asymmetries, we also studied the case with the event reconstruction turned on but 
with the trigger turned off. Similarly we also studied the case with the event reconstruction turned off 
but with the trigger turned on. The effect of energy smearing is also investigated. In each case, we 
computed the lepton correlated asymmetries, the neutrino correlated asymmetries and the W correlated 
asymmetries with respect to the planes described earlier. In Table |l|, we present the measured correlated 
asymmetries Ai,A 2 ,A 3 and A 4 for the four cases: (I), trigger off, reconstruction off; (II). trigger on, 
reconstruction off; (III), trigger off, reconstruction on; (IV). trigger on, reconstruction on; (V). trigger 
on, reconstruction on and with energy smearing included. The corresponding asymmetries when the top 
spins were "uncorrelated" (that is, their spins are summed over in their productions and averaged over in 
their decays,) are given in brackets. The average values and the standard deviations of the asymmetries 
were extracted directly from the simulated data. When both top and W spins were "uncorrelated", as 
in most standard event generator packages, we verified that all the asymmetries vanish if the trigger and 
event reconstruction were turned off. The top quark mass is taken to be m t — 176 GeV in all numerical 
analyses. As one can see in the table, the asymmetries measured by charged leptons are generally larger 
than the asymmetries measured by neutrinos and W bosons. The systematic effects of trigger and 
reconstruction are quite obvious in the case of neutrino asymmetries as one would expect from the large 
missing E±_ cut as well as the contributions from other neutrinos in the event. Due to space, a detailed 
discussions of the various effects and their origins will be given elsewhere. 

From the results of these simulation, we conclude that the lepton correlated asymmetries arising 
from the tt spin correlation can be observed at 3a level with a few hundred top dilepton events. This 
is certainly reachable with the projected luminosity of the Tevatron for Run II, with improved detector 
resolution and acceptance expected for both CDF and DO detectors upgrades and perhaps with improved 
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algorithms for identifying the top dilepton events. 

There are also plenty of rooms for improvement on our analysis of the correlated asymmetries. Even 
though trigger and reconstruction show little systematic effects for the lepton asymmetry, they do shift 
the neutrino and the W asymmetries by non-negligible ammounts. A quantitative understanding of 
these effects may allow us to tune the trigger selection and reconstruction algorithm for the observation 
of the asymmetries. It is well known that discrete ambiguities exist in general in the reconstruction of 
top dilepton events. A detailed quantitative study of its effect on the asymmetries could precipitate an 
improvement on reconstruction algorithm for better efficiency and better signal-to-noise ratio. We have 
been quite casual in choosing the planes to define the asymmteries and in choosing the combinations of 
these asymmetries to measure. One may wonder if there is an optimal choice of plane or combinations 
of planes that can maximize the observability of the correlation effect. One may also wonder if there are 
choices that will enhance or suppress the contribution of qq annihilation channel relative to the gg fusion 
channel. 

Before we conclude, we shall mention two other workfla which appear in the literature recently. These 
two papers ignore the off-diagonal terms of the top production density matrix and discuss only the spin 
correlated asymmetry. In contrary, our work addresses the more general density correlated asymmetry in 
top production. Only the contributions of the diagonal terms can be described as spin correlated. When 
off-diagonal terms are included in the analysis, most of the time the top quarks are produced in mix 
states of helicity(or spin) variable. As a result, our analysis includes effects that cannot be described as 
spin correlated asymmetry. Of course, for some observables the contributions of the off-diagonal density 
matrix may happen to be negligible. However it is generally not the case. In fact the optimal case we 
discovered, the Ai(l) asymmetry in the Table, cannot be described correctly by ignoring the off-diagonal 
terms. While the spin correlated asymmetry may enjoy simpler intuitive understanding, it, however, does 
not reflect the full range of possible interference effects emersed in the density matrix correlation in the 
top pair prodction. 

In conclusion, top spin correlation is certainly one of the most interesting top physics to be uncovered. 
It can provide a direct observation of the spin 1/2 character of top quark (which we have not been able 
to do for the lighter quarks) and can potentially test the V — A character of the weak charged current 
associated with top. We have clearly demonstrated the possibility of observing this effect at Tevatron 
Run II. The fact that this effect may be even harder for LHC to measure should make the task more 
important for Tevatron. A detailed account of our study will be presented elsewhertfl. 
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Table 1: The measured correlated asymmetries Aj(Z)'s, Ai(u)'s and Ai(W)'s for the four cases: (I), trigger off, reconstruc- 
tion off; (II). trigger on, reconstruction off; (III), trigger off, reconstruction on; (IV). trigger on, reconstruction on; (V). 
trigger on, reconstruction on with energy smearing. The corresponding asymmetries when the top spins were "uncorrelated" 
are given in brackets. All the values are in unit of percentage. The total number of simulated events is 100000. 
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